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Abstract 


Let Gr(d, n ) be the Grassmannian of d-dimensional linear subspaces of an ro-dimensional 
vector space V n . A submanifold A C Gr (d,n) gives rise to a differential system £(A) 
that governs d-dimensional submanifolds of V n whose Gaussian image is contained in X. 
Systems of the form £(A) appear in numerous applications in continuum mechanics, theory 
of integrable systems, general relativity and differential geometry, and include such well- 
known examples as the dispersionless Kadomtsev-Petviashvili equation, the Boyer-Finley 
equation, Plebansky’s heavenly equations, and so on. 

In this paper we concentrate on the particularly interesting case of this construction which 
however contains all essential challenges of the general theory, namely the case when A is a 
fourfold in Gr(3, 5). Our main goal is to investigate differential-geometric and integrability 
aspects of the corresponding systems £(A). There exist several approaches to dispersionless 
integrability: 

• The method of hydrodynamic reductions, based on a decomposition of a given multi¬ 
dimensional system into a collection of commuting (1 + l)-dimensional systems of hy¬ 
drodynamic type. Exact solutions defined by these reductions can be considered as 
natural dispersionless analogues of finite-gap solutions of integrable soliton equations. 

• The method of dispersionless Lax pairs, that is, pairs of Hamilton-Jacobi type equations 
whose compatibility conditions are equivalent to the system £(A). 

• Integrability on solutions , based on the requirement that the characteristic variety of 
an integrable system £(A) defines an integrable background geometry (Einstein-Weyl 
geometry in 3D, self-dual conformal geometry in 4D on every solution). 

• Integrability on equation , meaning integrability (in twistor-theoretic sense) of the canon¬ 
ical GL(2,M) structure induced on a fourfold X C Gr(3, 5). 

We demonstrate that for (d,n) = (3,5) these seemingly different approaches lead to one 
and the same class of integrable systems £(A) and prove that the parameter space of such 
systems is 30-dimensional. Factored by the natural action of the equivalence group SL(5), 
this gives 6-dimensional moduli space of integrable systems £(A). We give a complete de¬ 
scription of linearisable systems (the corresponding fourfold A is a linear section of Gr(3, 5)) 
and linearly degenerate systems (the corresponding fourfold A is the image of a quadratic 
map P 4 --■» Gr(3, 5)). The fourfolds corresponding to ‘generic’ integrable systems are not 
algebraic, and can be parametrised by generalised hypergeometric functions. 

MSC: 37K10, 37K25, 53A30, 53A40, 53B15, 53B25, 53B50, 53Z05. 

Keywords: Dispersionless Integrable System, Dispersionless Lax Pair, Hydrodynamic Re¬ 
duction, Einstein-Weyl Geometry, Submanifold of the Grassmannian, GL(2,R) Geometry. 
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1 Introduction 


1.1 Formulation of the problem 

In the most general setting the problem that we address in this paper can be described as fol¬ 
lows. Let Gr(d, n ) be the Grassmannian of d-dimensional linear subspaces of an n-dimensional 
vector space V n . A submanifold X C Gr(d, n) gives rise to the differential system E(X) that 
governs d-dimensional submanifolds of V n whose Gaussian image is contained in X (we recall 
that the Gaussian image of a submanifold is the collection of its tangent spaces translated to the 
origin). In this sense we are in the context of Grassmann geometries as discussed in [35]. Since 
d-dimensional submanifolds of V n are (locally) parametrised by n — d functions of d variables, 
we will assume that the codimension of X in Gr(d, n) also equals n — d: in this case E(X) will 
be a determined system of n — d first-order PDEs for n — d unknown functions of d independent 
variables. Systems of the form E(X) appear in a wide range of applications in differential geom¬ 
etry (in particular, calibrated geometries [35]), general relativity (Einstein-Weyl structures [18], 
heavenly-type equations [49]), and the theory of integrable systems (Backlund transformations 
in higher dimensions, dispersionless Lax pairs [61], dispersionless limits of multi-dimensional 
soliton equations). Particularly interesting examples are provided by multidimensional (d > 3) 
integrable systems E(X): such systems possess an infinity of multi-phase solutions that can be 
considered as analogues of multi-gap solutions of integrable soliton equations. It was observed 
in [23] that the requirement of the existence of such solutions constitutes an efficient criterion 
(known as the method of hydrodynamic reductions: see Sect. 1.5 for a brief description of the 
method), that allows one to classify integrable systems of the form £(X). In this paper we 
concentrate on the following key questions: 

• What is the dimension of the moduli space of integrable systems S(X)? What are the 
most interesting examples? 

• How does integrability of a system S(X) translate into differential geometry of the corre¬ 
sponding submanifold X C Gr(d, n)? 

In the case n = d + 1 the system E(X) consists of a single first-order PDE for a scalar 
function of d independent variables, which can be solved by the method of characteristics via 
reduction to ODEs. Thus, we can assume n > d+2. The case d = 1 also corresponds to an ODE 
system. The case d = 2, n arbitrary, leads to another familiar class of systems: introducing in 
V n coordinates t,x,u l , i = 1,..., n — 2, and parametrising two-dimensional submanifolds of V n 
in the form u l = u l (t,x), we can represent the corresponding system E(X) as 

F\ui,ui) = 0, ..., F- 2 K,<) = 0. 

Solving these equations for u\ in the form u\ = f l (u x ), here u = (u 1 ,... ,u n ~ 2 ), differentiating 
by x and setting u x = v, we obtain a system of conservation laws, 

= d x f{v). 

Differential-geometric and integrability aspects of such systems were extensively studied in [57, 
53, 17, 2], see also references therein. Thus, in what follows we assume n > d + 2 and d > 3. 

To be more specific we will concentrate on the particularly interesting case of this prob¬ 
lem where our results are fairly complete, the case of fourfolds X C Gr(3, 5) (n = 5, d = 3). 
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Introducing in V 5 coordinates x 1 , x 2 , x 3 , it, v, one can parametrise three-dimensional subman¬ 
ifolds of V 5 in the form u = u(x 1 ,x 2 ,x 3 ), v = v(x 1 ,x 2 ,x 3 ). Their tangent spaces are given 
by du = Uidx l , dv = Vidx * where Ui, Vi can be viewed as local coordinates on Gr(3,5). The 
corresponding system E(X) reduces to a pair of first-order PDEs for u and v, 

F(ui,U2,U 3 ,V 1 ,V 2 ,V 3 ) = 0 , G(ui,u 2 ,u 3 ,vi,v 2 ,v 3 ) = 0 , ( 1 ) 

Ui = du/dx 1 , Vi = dv/dx l . Equations (1) specify a fourfold X C Gr(3, 5). The class of systems 
( 1 ) is invariant under the equivalence group SL(5) that acts by linear transformations on the 
combined set of variables x 1 , x 2 , x 3 , u, v. Since this action preserves the integrability, all our 
classification results will be formulated modulo SL(5)-equivalence. Necessary details on the 
equivalence group are provided in Sect. 3. 

There exists a whole variety of important examples that fall into class (1). One of them ap¬ 
pears in the context of the dispersionless Kadomtsev-Petviashvili (dKP) equation, u x t — u x u xx — 
Uyy = 0, one of the most well-studied dispersionless integrable PDEs arising in nonlinear acous¬ 
tics [59] and the theory of Einstein-Weyl structures [18]. Its dispersionless Lax pair [61] consists 
of two first-order relations of type (1), 

1 2 1 .. 

v y - -v x -u x = 0, v t - -v x - v x u x -u y = 0, (2) 

here (a: 1 , x 2 , x 3 ) = ( x,y,t ). The dKP equation results from (2) on elimination of v, that is, via 
the compatibility condition v y t = vt y . Similarly, the elimination of u leads to the modified dKP 
(mdKP) equation, v x t — (v y — \v 2 )v xx — v yy = 0. Thus, relations (2) provide Backlund-type 
transformation connecting dKP and mdKP equations. We refer to Sect. 4 for further examples 
and classification results. 


1.2 Non-degeneracy condition 

In what follows we assume that system (1) is non-degenerate in the following sense. Let us 
consider the corresponding linearised system, 

O £■)(“) +o £-)(") +(£> J0(v) 

V Lr Vl J \ V J x \ \ ' J ’V 2 J \ V J x 2 V ^ u 3 ^^3 J V * J x 3 

obtained by setting u —> u + eW, v —> v + eV, expanding F and G in Taylor series and keeping 
terms of the order e. Non-degeneracy means that the dispersion relation (characteristic variety), 


det 



F u i Fvi \ i \2 f F U2 F V2 
G U1 G V1 J ' A V G U2 G V2 


+ A 3 


F U3 F V3 
Gu 3 G V3 


= 0 , 


defines an irreducible conic in IP 2 with homogeneous coordinates (A 1 : A 2 : A 3 ). Explicitly, the 
dispersion relation can be represented in the form (A 1 , A 2 , A 3 )< 7 **(A 1 , A 2 , A 3 ) 4 = 0 where is the 
3x3 symmetric matrix, 

= gV = ^ (F u i G v j + F u j G v i - F v i G u j - F v j G u i ). 

It gives rise to the conformal structure g = gijdx l dx J (here gij is the inverse of g lJ ). Note 
that non-degeneracy is equivalent to det < 7 ^ 0 : this is the case for all known systems of physi¬ 
cal/geometric relevance. It turns out that the signature of g is always Lorentzian, and thus our 
PDE system is hyperbolic. 
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Geometric aspects of conformal structures defined by the characteristic variety will play a 
key role in our characterisation of integrable systems: we will see that solutions to integrable 
equations carry ‘integrable’ background geometry. In 3D, this is the Einstein-Weyl geometry. 

1.3 Einstein-Weyl geometry in 3D 

Recall that an Einstein-Weyl structure consists of a symmetric connection B and a conformal 
structure g such that: 

(a) connection B preserves the conformal class of g: B[g] = 0, 

(b) trace-free part of the symmetrized Ricci tensor of B vanishes. 

In coordinates, this gives 


^k9ij — ^kdiji R(ij) — A.9ij i (3) 

where u = u>kdx k is a covector, R(ij) is the symmetrized Ricci tensor of B, and A is some function 
[12], Note that it is sufficient to specify g and u only, then the first set of equations uniquely 
determines B. The integrability of Einstein-Weyl equations (3) by twistor-theoretic methods 
was established by Hitchin [36]. It was shown in [22] that generic Einstein-Weyl structures 
are governed by the Manakov-Santini system introduced in [39] as a two-component integrable 
generalisation of the dKP equation. We will see that solutions to integrable systems (1) carry 
Einstein-Weyl geometry: the conformal structure g defined by the characteristic variety must 
be Einstein-Weyl on every solution, furthermore, the covector u: can be expressed in terms of g 
by the universal formula 

Wfc = Zgkj’Dx* ( 9 js ) + T> x k (In det gij), (4) 

where T> x k denotes total derivative with respect to the independent variable x k (note that g 
depends on first-order jets of the solution u, v). Formula (4) appeared in [27] in geometric 
approach to the dispersionless integrability in 3D. It is invariant under the gauge transformation 
g —>• A g, uj —> oj + din A, the property characteristic of Einstein-Weyl geometry. According to 
the result of Cartan [12], the Einstein-Weyl property of a triple (B, g, u) is equivalent to 
the existence of a two-parameter family of surfaces that are null with respect to the conformal 
structure g (that is, tangential to the null cones of g), and totally geodesic in the Weyl connection 
B. In the context of integrable systems (1), such surfaces are provided by the corresponding 
dispersionless Lax pairs: these consist of A-dependent vector fields X, Y that are required to 
commute modulo (1), identically in the ‘spectral parameter’ A. For systems (1), the existence of 
such Lax pairs is equivalent to the integrability by the method of hydrodynamic reductions, see 
Sect. 2.2. Taking integral surfaces of the distribution spanned by X, Y in the extended four-space 
with coordinates x , y, t, A, and projecting them down to the space of independent variables x, y , t, 
we obtain the required two-parameter family of null totally geodesic surfaces. Let us mention 
that relations of dispersionless integrable systems in 3D to Einstein-Weyl geometry have been 
discussed previously in [58, 11, 19, 27], see also references therein. 

Example 1. Let us consider the system 

1 9 

U t - -u x - Vy = 0, v x -Uy = 0, 

which reduces to the dKP equation, u x t — u x u xx — u yy = 0, on elimination on v. Its characteristic 
variety defines the conformal structure g = 4dxdt — dy 2 + 4 u x dt 2 . Introducing the covector ui = 
—Au xx dt by formula (4), one can verify that the pair g,u satisfies the Einstein-Weyl equations if 
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u solves the dKP equation. This Einstein-Weyl structure was obtained previously in [18]. The 
corresponding Lax pair has the form 

A — dy T Uxx&X: ^ (A T "P i^xx^ "P ^xy)^ A? 

one can verify that these vector fields commute modulo the above system. Projecting integral 
surfaces of the distribution spanned by X, Y from the extended space of variables x , y , t, A to the 
space of independent variables x , y , t, one obtains a two-parameter family of null totally geodesic 
surfaces of the corresponding Einstein-Weyl structure. 

1.4 GL(2,M) structures 

The tangent bundle to the Grassmannian Gr(3, 5) carries canonical generalised conformal struc¬ 
ture defined by the family of Segre cones duidvj — dujdvt = 0. Thus, each projectivised tangent 
space PTGr(3,5) contains a Segre variety, an algebraic threefold of degree three. Given a non¬ 
degenerate fourfold X C Gr(3, 5) (a fourfold is said to be non-degenerate if it gives rise to a 
non-degenerate system (1)), the intersection of its tangent space TA - with the Segre cone is a 
two-dimensional rational cone of degree three; its projectivisation is a rational normal curve of 
degree three (twisted cubic). Thus, PTAT is supplied with a field of twisted cubics. This is known 
as a GL(2,M) structure on X. It was demonstrated by Bryant [8] that every four-dimensional 
GL(2,M) structure defines on X a canonical affine connection which preserves the GL(2,M) 
structure, and whose torsion lies in 8-dimensional irreducible representation of GL(2,R). We 
will call it the Bryant connection, see Sect. 5.4 for computational formulae. Various important 
properties of system (1) have natural interpretation in terms of this connection: 

• System (1) is linearly degenerate and integrable if and only if the Bryant connection is 
symmetric and flat (Proposition 5 of Sect. 4.2). 

• System (1) is integrable if and only if the curvature R and the covariant derivative VT of 
the torsion T of the Bryant connection are certain invariant quadratic expressions in T, 

R = f(T 2 ), XT = g(T 2 ), 

see Theorem 3 of Sect. 5.4 for precise statements. These expressions are analogous to the 
ones obtained by Smith [55] in the context of 5-dimensional GL(2, M) structures associated 
with integrable equations of the dispersioless Hirota type [26]. Our expressions provide a 
compact invariant formulation of the integrability conditions. 

1.5 The method of hydrodynamic reductions 

In the most general set-up, the method of hydrodynamic reductions [23] applies to quasilinear 
systems of the form 

d(u)u x + B(u)u y + C(u)u t = 0, (5) 

where u = (u 1 , ...,u m ) t is an m-component column vector of the dependent variables, and 
A,B,C are l x m matrices where Z, the number of equations, is allowed to exceed the number 
of unknowns, m. Note that system (1) can be brought to form (5) by choosing Ui,Vi as the new 
dependent variables and writing out all possible consistency conditions among them - see Sect. 
2.1 for details. The method of hydrodynamic reductions consists of seeking multiphase solutions 
in the form 

u = u (R 1 ,..., R n ), 
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where the phases R l (x,y,t), whose number N is allowed to be arbitrary, are required to satisfy 
a pair of consistent (1 + l)-dimensional systems, 

K = A \R)K, R'y = fj(R)R* x , (6) 

known as systems of hydrodynamic type. The corresponding characteristic speeds A* and /T are 
required to satisfy the commutativity conditions [57], 

djA* = dj/S 

\i - A* fil - n l ’ 1 ' 

here i ^ j, dj = d^j. Solutions of this type originate from gas dynamics, and are known as 
nonlinear interactions of planar simple waves. They can be interpreted as natural dispersionless 
analogues of finite-gap solutions of 2 + 1 dimensional soliton equations. Equations (6) are said 
to define an IV-component hydrodynamic reduction of the original system (5). System (5) is 
said to be integrable if, for every N, it possesses infinitely many IV-component hydrodynamic 
reductions parametrised by N arbitrary functions of one variable. This requirement imposes 
strong constraints (integrability conditions) on the matrix elements of A, B and C. We refer to 
[23] for further details and references. 

1.6 Summary of the main results 

In Sect. 2.1 we apply the method of hydrodynamic reductions to equations of type (1). This 
results in an overdetermined involutive system of PDEs for the functions F and G (the integra¬ 
bility conditions). The analysis of this system leads to our first result (see Theorem 1 of section 
2.1 and Proposition 2 of section 3): 

The moduli space of non-degenerate integrable systems (1) is 6-dimensional. Furthermore, 
the following conditions are equivalent: 

(a) System (1) is integrable by the method of hydrodynamic reductions. 

(b) Conformal structure g defined by the characteristic variety of system (1) is Einstein-Weyl 
on every solution (with covector u given by formula (4))- 

Thus, Einstein-Weyl equations (3) provide an efficient integrability test. Using the integra¬ 
bility conditions we prove that every integrable system (1) possesses a dispersionless Lax pair 
(Theorem 2 of Sect. 2.2). Furthermore, we describe a construction that links Lax pairs to null 
totally geodesic surfaces of the corresponding Einstein-Weyl structure (Sect. 2.3). 

In Sect. 3 we demonstrate that the class of equations (1) is invariant under the equivalence 
group SL(5) that acts by linear transformations on the combined set of variables x 1 , x 2 , x 3 , u, v. 
All our results are formulated modulo this equivalence. 

In Sect. 4 we describe a variety of integrable examples of type (1) expressible in elementary 
functions, modular forms and theta functions. In particular, we demonstrate that system (1) is 
linearisable by a transformation from the equivalence group if and only if the conformal structure 
defined by its characteristic variety is flat on every solution (Sect. 4.1). Next, we show that 
linearly degenerate integrable systems are characterised by the property that the associated 
Bryant connection is symmetric and flat. Fourfolds X corresponding to such systems are images 
of quadratic maps P 4 —-> Gr(3, 5) given by the classical construction of Chasles. This provides 
a complete list of normal forms of linearly degenerate integrable systems (Sect. 4.2). Based on 
the construction of Odesskii and Sokolov [44, 45] we point out that generic integrable systems 
(1) can be parametrised by generalised hypergeometric functions (Sect. 4.4). This demonstrates 




that fourfolds X C Gr(3, 5) corresponding to integrable systems (1) may have a fairly intricate 
analytic structure: in particular, they do not need to be algebraic. Among the simplest nontrivial 
examples of this kind one can mention the system 

Uy \ . \ 9 

Ut = V X , v t = — + -r]{u x )v x , 
v x 6 

which appeared in [46] in the classification of integrable hydrodynamic chains. In this case the 
integrability conditions result in the Chazy equation for 77 , rj'" + 2 rjrj" — 3r]' 2 = 0 (here prime 
denotes differentiation by u x ), whose generic solution is expressible in terms of the Eisenstein 
series E^. 

In Sect. 5 we demonstrate that every fourfold X C Gr(3, 5) inherits a GL(2,R) struc¬ 
ture, namely, a field of twisted cubics specified in the projectivised tangent bundle PTX (Sect. 
5.1). The integrability can be reformulated geometrically as the requirement that the associ¬ 
ated GL(2,R) structure possesses infinitely many three-dimensional submanifolds carrying a 
holonomic characteristic net (Sect. 5.2). We show that every fourfold X corresponding to in¬ 
tegrable system (1) carries canonical conformally symplectic structure which is parallel in the 
Bryant connection (Sect. 5.3). Finally, we reformulate the integrability conditions of system (1) 
in terms of the curvature and torsion of the associated connection (Theorem 3 of Sect. 5.4). 
We emphasize that GL(2,R) structures arising in this context are quite different from those 
discussed in [8], in particular, they generally have nonzero torsion. 

Throughout the paper, most of our considerations are local, although some results related to 
algebraic fourfolds X C Gr(3,5) have global nature. We are primarily interested in real smooth 
structures and hyperbolic systems, however, the theory applies to the complex-analytic case 
as well. In the calculations we use computer algebra systems Mathematica and Maple (these 
only utilise symbolic polynomial algebra over Q, so the results are rigorous). The programs are 
available from arXiv:1503.02274vl [math.DG]. 

2 Integrability conditions and Lax pairs 

2.1 Derivation of integrability conditions 

In this section we establish our first main result: 

Theorem 1 The parameter space of non-degenerate integrable systems (1) is 30- dimensional. 
Furthermore, the following conditions are equivalent: 

(a) System (1) is integrable by the method of hydrodynamic reductions. 

(b) Conformal structure g defined by the characteristic variety of system (1) is Einstein-Weyl 
on every solution (with covector 10 given by formula (4))- 

Proof: 

Applied to system (1), the method of hydrodynamic reductions leads to a set of differen¬ 
tial constraints (integrability conditions) for the functions F and G. The 30-dimensionality of 
the parameter space follows from the involutivity of these conditions. Once the integrability 
conditions are derived, the equivalence of (a) and (b) can be shown by a direct calculation. 
Alternatively, one can utilise the fact that every integrable system (1) possesses a dispersionless 
Lax pair (Theorem 2 of Sect. 2.2), which provides null totally geodesic surfaces of the corre¬ 
sponding Einstein-Weyl structure (O, g,w). This implies the Einstein-Weyl property due to the 
result of Cartan [12]. 
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To derive the integrability conditions we proceed as follows. Let us first rewrite system (1) 
in evolutionary form, 


u t = f{u x ,Uy,v x ,v y ), v t = g{u x ,Uy,v x ,v y ), (8) 

here we use the notation (x 1 ,^ 2 ,^ 3 ) = ( x,y,t ). Our strategy is to derive a set of constraints 
for the right hand sides / and g that are necessary and sufficient for integrability. We begin 
by transforming our system into first-order quasilinear form. This can be done by taking first- 
order partial derivatives of u and v as the new dependent variables and writing out all possible 
consistency conditions among them. Thus, we introduce the notation 


u x = a , u y = b , v x =p, v y = q, u t = f(a,b,p,q), v t = g(a,b,p,q). 


This provides an equivalent quasilinear representation, 

a y = b x , a t = f (a, b,p, q) x , b t = f{a,b,p,q) y , 

(9) 

Py = Qx , Pt = g(a,b,p,q) x , q t = g{a,b,p,q) y . 

Note that quasilinear representation (9) is of type (5) with m = 4, 1 = 6, u = ( a,b,p,q ). 
Looking for multi-phase solutions in the form 

a = a(R 1 ,...,R N ), b = b(R x , ...,R N ), p = p{R l , ...,R N ), q = q{R l ,..., R N ), 

where the phases R l satisfy equations (6), and substituting this ansatz into (9), we obtain the 


relations 



dib = p l did, 

d%q = pfdip, 

(10) 

here di = d R i, as well as 



(A* -fa~ P l h)did = {f p + p l fq)dip, 

(A* - g P - p l g q )dip = {g a + g l gb)did. 

(11) 


The last two equations imply the dispersion relation connecting A* and //, 

(A* -fa- g l fb)(X l -g P - g l g q ) = {f P + g l f q )(g a + g l gb)- ( 12 ) 


In what follows we assume that the dispersion relation defines an irreducible conic in the (A, p)- 
plane: this is equivalent to the non-degeneracy of system (1) as discussed in Sect. 1.2. Setting in 
(11) did = (p l dip and solving the resulting equations for A* and p l we can parametrise dispersion 
relation (12) in the form 

i = fp + (fa ~ 9p V ~ 9a y = (fq + /fe^)(gp + 9a^) ~ (fp + fa^iOg + Qb^) 

fq + (fb - gqW - gbP i2 ’ fq + (fb ~ W ~ 9b^ 

Substituting these parametric expressions into the commutativity conditions (7), and using the 
relations 

did = p l dip, dib = p l ip l dip, dtq = g l d t p, (13) 

we obtain djp 1 in the form djp 1 = (... )djp, i j, where dots denote rational expressions 

in < p l , <p J whose coefficients depend on second-order partial derivatives of / and g (we omit 

these expressions due to their complexity). Calculating consistency conditions for relations 


10 




(13), didja = djdid, didjb = djdib, didjq = djdiq, we obtain (one and the same!) expression for 
didjp in the form didjp = (... )dipdjp, i / j, where, again, dots denote terms rational in p> 1 and 
ip 1 . Ultimately, IV-phase solutions are governed by the relations 

dj<p l = (■■■ )djp, didjp = (... )dipdjp, i + j. (14) 

To ensure their solvability we need to impose the compatibility conditions 

dkdjip 1 = djd k ip\ d k djdip = djd k dip, (15) 

that are required to hold for every triple of indices i ^ j ^ k. Direct calculation based on (13) 
and (14) results in 

dkdjtp 1 - djd k <p l = (... )djpd k p , d k djdip - djd k dip = (... )d i pd j pd k p 1 

where dots denote rational expressions in tp l , p 3 , ip k whose coefficients depend on partial deriva¬ 
tives of / and g up to order three. To ensure solvability of equations (14) we set all these 
coefficients equal to zero. This gives a system of differential constraints for / and g that is linear 
in the third-order derivatives thereof. Moreover, the number of linearly independent compati¬ 
bility conditions equals the total number of third-order partial derivatives of / and g, namely 
40. Solving for third-order derivatives of / and g we obtain the required integrabililty conditions 
that can be represented in symbolic form, 

d 3 f = R(df , dg, d 2 f , d 2 g), d 3 g = S(df, dg, d 2 f, d 2 g), (16) 


40 equations altogether. Here R and S depend rationally on the first- and second-order partial 
derivatives of / and g. Straightforward calculation shows that overdetermined system (16) is 
in involution. Thus, the parameter space of integrable systems (1) is 30-dimensional: one can 
arbitrarily specify the values of f,g , df,dg and d 2 f,d 2 g at any fixed point. This amounts to 
2x1 + 2x4 + 2x10 = 30 arbitrary constants. 

Finally, the equivalence of (a) and (b) can be established by a direct calculation: in coordi¬ 
nates ( x,y,t ), the characteristic variety of system (8) defines contravariant metric (metric with 
upper indices) with the matrix 


g lJ = 


( fa9p fp9a 2 if a 9q fq9a T fb9p fp9b ) 

2 (fa,9q fq9a + fb9p ~ fp9b ) fb9q ~ fq9b 


2 if a 4 “ 9p) 


-\{fb+9q) 


~\{fa + 9p) \ 

~\{fb + 9q) 

1 / 


Introducing uj by formula (4) and calculating Einstein-Weyl conditions (3) we obtain expressions 
depending on third-order partial derivatives of u and v. These have to vanish identically modulo 
(8), in other words, on every solution. Using (8) and its differential consequences to eliminate 
all partial derivatives involving differentiation by t, we obtain differential expressions that are 
polynomial in the remaining second- and third-order partial derivatives of u and v. It can 
be verified directly that the vanishing of all coefficients of these polynomials is equivalent to 
integrability conditions (16). More precisely, the Einstein-Weyl conditions contain terms of 
two types: linear in third-order derivatives, and quadratic in second-order derivatives of u and 
v. Coefficients at third-order derivatives of u and v vanish identically due to the choice of oj, 
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while coefficients at quadratic terms (105 coefficients altogether) give all of the 40 integrability 
conditions (16). Details of calculations described above, including the integrability conditions, 
are available from arXiv:1503.02274vl [math.DG]. This finishes the proof of Theorem 1. 

Remark. Note that each compatibility condition (15) involves three distinct indices only. This 
leads to the following important observations. 

(i) Every system of type (1) has infinitely many two-component reductions parametrised by 
two arbitrary functions of one variable. Indeed, these reductions are governed by equations (14) 
where i,j = 1,2. Since in this case equations (14) are automatically consistent, their general 
solution depends, modulo reparametrisations R l tp l (R l ), on two arbitrary functions of one 
variable. Therefore, the existence of two-component reductions is a common phenomenon that 
is not related to the integrability. 

(ii) On the contrary, the existence of three-component reductions implies the existence of N- 
component reductions for arbitrary N. Thus, one can define integrability as the existence of 
infinitely many three-component reductions parametrised by three arbitrary functions of one 
variable. This property is reminiscent of the well-known three-soliton condition in the theory of 
(1 + l)-dimensional integrable soliton equations. 

2.2 Dispersionless Lax pairs 

We say that system (1) possesses a dispersionless Lax pair if there exists an auxiliary function 
S satisfying two Hamilton-Jacobi type equations, 

S 2 = P(Si, Ui,Vi), S 3 = Q(Si,Ui,Vi), (17) 

such that the compatibility condition, S 23 = £ 32 , holds identically modulo (1); here S) = dS/dx 1 . 
Note that the dependence of P and Q on S 1 is usually nonlinear. Lax pairs of this type appeared 
in the context of dispersionless integrability in [61]; in many cases, including the one studied 
in this paper, the existence of a dispersionless Lax pair is equivalent to the integrability by the 
method of hydrodynamic reductions. For instance, system (2) possesses the Lax pair 

S y = S% + v x S x , St = —S x + 2 'v x S x + ( u x + v x )S x ; 

the compatibility condition S y t = St y is satisfied identically modulo (2). The main result of this 
section is the following: 

Theorem 2 Integrability of system (1) is equivalent to the existence of a ‘generic’ dispersionless 
Lax pair. 

The definition of the term ‘generic’ is given after formula (19). 

Proof: 

Let us again consider system (1) in evolutionary form, 

— f j 'Uy j Lr j Vy ) > j Hy j Lr j Vy ) ■ 

Setting u x = a, u y = b, v x = p, v y = q we will look for a Lax pair in the form 

S y = P(S x ,a,b,p,q), S t = Q{S x ,a,b,p,q). (18) 
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Requiring that the consistency condition S y t = Sty holds identically modulo (9), we obtain 6 
first-order relations for P and Q: 


faPa + 9aPp P P\Qa ~ Q\Pa = 0, fpP a + 9pPp + P\Qp ~ Q\Pp = 0, 
Qa = fbPa + faPb + 9bPp P 9a,Pq P P\Qb ~ Q\Pb, 

Qp = fqPa P fpPb P 9qPp “1“ 9pPq P P\Qq Q\Pqi 


(19) 


Qb — fbPb p 9bPqi Qq — fqPb “I - 9qPqi 

here A = S x . We say that system (1) possesses a ‘generic’ Lax pair if relations (19) are in 
involution, and P,Q do not satisfy any extra first-order relations other than (19). Generic Lax 
pair is a generic solution to (19). Differentiating each of these relations by a,b,p,q and A, we 
obtain 30 relations which, in the non-degenerate case, can be uniquely resolved for all second- 
order partial derivatives of P and Q, thus leading to a closed system. It can be verified directly 
that the resulting system is involutive if and only if the functions / and g satisfy integrability 
conditions (16). This finishes the proof of Theorem 2. 

Remark 1. The general solution of the involutive system for P and Q depends on 6 arbitrary 
constants. This comes from the count 2 + 10 — 6 = 6, where 2 corresponds to the values of P and 
Q, 10 comes from the values of their first-order derivatives, and —6 is due to the 6 first-order 
relations (19). On the other hand, the Lax pair is form-invariant under transformations 

S aS + (3u P 'yv + g,x + uy + rjt, 


which also involve 6 arbitrary constants. Modulo these transformations, ‘generic’ Lax pair is 
essentially unique. 

Remark 2. One can show that, as a consequence of relations (19), the quantities P\ and Q\ 
satisfy the dispersion relation, 


det 


fa fp 

9a 9p 


p 


fb fq 

9b 9q 


P\~ 


1 0 
0 1 


Q\ 


= 0. 


( 20 ) 


This relation plays an important role in the geometric interpretation of Lax pairs discussed 
below. 


2.3 Lax pairs and null totally geodesic surfaces 

Here we outline a general construction that leads from Lax pair (18) to null totally geodesic 
surfaces of the corresponding Einstein-Weyl structure (ID, g, uj). Differentiating (18) by x and 
setting S x = A we obtain 

A y — P\K P Pa&x P Pb^x P PpPx P PqQxi A t — Q\^x T Qa&x T Qb^x P QpPx P QqQx• (21) 

With this system we associate two vector fields, 

-^ = — P\~!b P {PaCLx p Pbbx p PpPx + Pq9x)-§\1 

A P A a 

Y — (Qa^x Qb^x QpPx QqQx)"Q\i 
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which live in the extended four-dimensional space with coordinates x, y, t, A. Note that the 
compatibility condition, A y t = A t y , is equivalent to the commutativity of these vector fields: 
[X, Y] = 0. The geometry behind this construction is as follows. Let us consider the cotangent 
bundle Z of a solution u(x,y,t ), v(x,y,t), with local coordinates (x,y,t,S x ,S y ,St). Equations 
(18) specify a four-dimensional submanifold M 4 C Z parametrised by x,y,t and A. The com¬ 
patibility of equations (18) indicates that this submanifold is coisotropic. Vector fields X, Y 
generate the kernel of the restriction to M 4 of the symplectic form dS x A dx + dS y A dy + dSt A dt. 
Equations (21) mean that X,Y are tangential to the Lagrangian submanifold of M 4 defined 
by the equation A = A (x,y,t). Solutions S(x,y,t ) of equations (18) can be interpreted as gen¬ 
erating functions of Lagrangian submanifolds contained in M 4 . Such submanifolds depend on 
one arbitrary function of a single variable; they can be obtained by taking any one-parameter 
subfamily of the two-parameter family of integral surfaces of the distribution (X, Y) foliating 
M 4 . 

Projecting the two-parameter family of integral surfaces of the distribution (X, Y) from A/ 4 
to the space of independent variables x,y,t we obtain a two-parameter family of null totally 
geodesic surfaces of the Weyl connection B. Indeed, projecting X and Y gives vector fields 

x - 9 p 9 v - 9 n 9 

ay ox at ox 

which commute if and only if A satisfies equations (21). It remains to show that (X,Y) is a 
null distribution (that is, tangential to the null cones of g ), and that the covariant derivatives 
B^X, B^y, ByX, ByY belong to ( X,Y). Equivalently, one can introduce the covector 8 = 
dx+P\dy+Q\dt that annihilates X, Y , and verify that 8 is null, and that B>^8A8 = B yd/\8 = 0. 
This follows from equations (19) satisfied by the functions P(A, a,b, c) and Q(X,a,b,c). In 
particular, the fact that 8 is null follows from identity (20). 

3 SL(5) as the equivalence group: geometry of Gr(3, 5) 

Let V 5 be 5-dinrensional vector space with coordinates x 1 , x 2 , x 3 , u, v, set p = (u,vY and x = 
(x 1 , x 2 , x 3 y. Three-dimensional linear subspaces of V are defined by the equation dp = t/dx 
where 

U = ( u i ^2 u 3 \ 

V vi v 2 v 3 ) ' 

Thus, the Grassmannian Gr(3,5) is six-dimensional, with an affine chart identified with the 
space of 2 x 3 matrices U. Solutions to system (1), ^(x 1 ,^ 2 ,® 3 ) and u(x 4 ,x 2 ,x 3 ), can be 
interpreted as three-dimensional submanifolds of V 5 . Their tangent spaces are specified by 
matrices U whose entries coincide with first-order derivatives of u and v: Ui = du/dx ®, Vi = 
dvjdx 1 . Thus, system (1) defines a fourfold X C Gr(3,5), while solutions to (1) correspond to 
three-dimensional submanifolds of V 5 whose Gaussian images are contained in X. The action 
of SL(5), 

( dp \ (A B\( dp A 

V dx ) \C D ) V dx ) ’ 

naturally extends to Gr(3, 5): 

U = (AU + B)(CU + D)- 1 ; (22) 
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here A, B , C, D are 2x2, 2x3, 3x2 and 3x3 matrices, respectively; notice that the extended 
action is no longer linear. Transformation law (22) suggests that the action of SL(5) preserves 
the class of equations (1), indeed, first-order derivatives transform through first-order derivatives 
only. Transformations of this form preserve integrability, so that the group SL(5) can be viewed 
as the natural equivalence group of the problem: it maps integrable systems to integrable systems. 
Thus, SL(5) is the point symmetry group of integrability conditions (16) derived in Sect. 2.1. 
The classification of integrable systems (1) will be performed modulo this equivalence: two 
SL(5)-related equations should be regarded as ‘the same’. 

Geometrically, our problem is reduced to the classification of fourfolds X of the Grassmannian 
Gr(3, 5) that satisfy certain ‘integrability’ conditions (to be specified later, see Sect. 5.2), modulo 
the action of SL(5). 

Given two infinitesimally close three-dimensional linear subspaces of V 5 defined by 2 x 3 
matrices U and U + dU, the condition that their intersection is two-dimensional is given by 
rank dU = 1. This condition is invariant under the action of SL(5) as specified by (22). Thus, 
each tangent space of Gr(3,5) is equipped with the cone C defined by the equation rank dU = 1, 
that is, by the vanishing of all 2 x 2 minors of dU. Projectivisation of this cone is known as 
the Segre variety; it is a non-singular algebraic threefold of degree three. The field of Segre 
cones supplies Gr(3, 5) with the generalised flat conformal structure; it is manifestly invariant 
under the action of SL(5). The converse is also true: transformations from SL(5) are the only 
diffeomorphisms that preserve the field of Segre cones. 

Proposition 1. The group of conformal automorphisms of the field of Segre cones is isomorphic 
to SL(5). 

This is a well-known fact: direct proof would consist of the calculation of conformal automor¬ 
phisms of the family of Segre cones. Let us point out that, in coordinates Ui, V{, the infinitesimal 
generators corresponding to equivalence transformations (22) are as follows: 

6 translations: 

U i = -S~, V * 

OUi OVi 


12 linear transformations (note the relation X„; = Lu + B 22 ): 


yr - 9 I d _ d d 

-'Mjj — Q H" i -L'll — q_. 5 -^12 — ^k 


s A 1 q 5 —-L-L '"K r\ 

OUj OVj OUfc OVh 

6 non-linear (projective) transformations: 


d d 

L 21 — Vk ~—, L 2 2 — vj .——. 

OUk OVk 


M d d d d 

Pi = UiUk 7T-1- ViUk , Q i = UiVk W-h ViVk • 

OUk OVk OUk OVk 


The ideal 12 defining the field of Segre cones is generated by quadratic forms, 

= span {duidvj — dviduj}, 

which are nothing but second fundamental forms of the Pliicker embedding of the Grassmannian, 
see e.g. [34]. It remains to point out that any vector held X satisfying Lxkl = 0 (mod 12) is 
spanned by the above infinitesimal generators. We refer to [33, 4, 31, 32] for generalisations of 
this Liouville-type result. 
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Consider the action of the equivalence group SL(5) on the space J 1 (M 4 , R 2 ) of 1-jets of func¬ 
tions /, g of variables a, b,p, q. This is a 14-dimensional space with coordinates Ui,Vi , /, g , f Ui , f Vi , 
g Ui ,9vi , i = 1, 2, which can be viewed as an affine chart in the bundle of 4-dimensional tangent 
subspaces of Gr(3, 5). The action of SL(5) is canonically defined in the latter space, but by 
abuse of notation we will be working with J 1 . 

Lemma. The group SL(5) has a unique Zariski open orbit in the space J 1 (R 4 ,R 2 ) (its comple¬ 
ment consists of 1-jets of degenerate systems). 

The proof of this statement is as follows. The group SL(5) acts transitively on the Grass- 
mannian Gr(3, 5) with the stabilizer of a point o being the parabolic subgroup P Q = S(GL(2) x 
GL(3)) x (R 2 (g)M 3 ) of upper-triangular block matrices of the size 2+3. As only S'(GL(2) x GL(3)) 
acts on T 0 Gr(3,5), (i.e. R 2 <S> R 3 acts trivially), the action is transitive on 4-planes correspond¬ 
ing to non-degenerate 1-jets of f,g. At the level of Lie algebra 51(5), prolongation of the above 
infinitesimal generators to J 1 (R , R 2 ) has full rank in the Zariski open set of non-degenerate 
1-jets. Indeed, the 24 x 14 matrix of coefficients of these vector fields drops rank precisely on 
the submanifold det g lJ = 0 where g 13 is the conformal structure defined by the characteristic 
variety. 

Remark 1. The above lemma allows one to assume that all sporadic factors depending on 
first-order derivatives of / and g that arise in the process of Gaussian elimination in the proofs 
of our main results, are nonzero. This considerably simplifies the arguments by eliminating 
unessential branching. Furthermore, in the verification of polynomial identities involving first 
and second-order partial derivatives of / and g one can, without any loss of generality, give the 
first-order derivatives any ‘generic’ numerical values: this often renders otherwise impossible 
computations manageable. 

In fact, even the action of SL(3) C P 0 on T 0 Gr(3, 5) is already transitive, and the stabilizer of 
a non-degenerate four-plane in T 0 Gr(3, 5) is the subgroup H a = GL(2) x (R 2 <8>R 3 ) of P 0 , where 
we identify R 3 ~ 5 2 (R 2 ) as an SL(2)-representation. Prolongation of this action to J 2 (R , R 2 ) 
is locally free near a generic point (straightforward calculation shows that the prolongation of 
infinitesimal generators of sl(5) to the space of 2-jets has full rank 24 at generic points). 

Proposition 2. The action of the equivalence group SL(5) on the 30 -dimensional parameter 
space of integrable systems (1) is locally free (in the Zariski open set of generic points). 

That is, ‘generic’ orbits of this action are 24-dimensional, so that ‘generic’ integrable systems 
of type (1) do not possess any continuous point symmetries coming from the equivalence group. 
Thus, the moduli space of integrable systems depends on 30 — 24 = 6 essential parameters. To 
prove the above statement let us note that the submanifold in the space of 3-jets of /, g given 
by integrability conditions (16) projects diffeomorphically onto the space of 2-jets J 2 (R 4 ,R 2 ) 
(which has dimension 30). 

Remark 2. The action of SL(5) on the 30-dimensional parameter space M 30 is algebraic, so 
that by the Rosenlicht theorem [51] there exists a quotient M r 3 g/SL(5), which is a rational 
algebraic variety of dimension 6 (the geometric quotient M 30 /SL(5) is singular). Indeed, the 
equivalence group transforms the initial conditions of system (19) rationally. 

4 Examples and classification results 

This section contains various classification results based on integrability conditions (16). 


16 


We produce an abundance of non-trivial examples of integrable systems of type (1), both known 
and new, expressible in elementary functions, theta functions, modular forms and generalised 
hypergeometric functions. 

4.1 Linearisable systems 

In this section we characterise systems (1) that can be linearised by a transformation from 
the equivalence group SL(5). Taking a linear system, say u 3 = v\, v 3 = u 2 , and applying 
transformations from the equivalence group, one obtains systems of Monge-Ampere type, 

a lj {uiVj - UjVi ) + b l Ui + c l Vi + m = 0 , , . 

a lj (UiVj — UjVi) + /3 l Ui + 7 l Vi + /i = 0 , 

where each equation represents a linear combination of minors of the 2x3 matrix 

JJ = ( “1 U 2 U 3 \ 

v Vi v 2 V 3 ) ' 

Conversely, in three dimensions every system of this form is linearisable. 

Proposition 3. For a non-degenerate system of type (1), the following conditions are equivalent: 

(a) System is linearisable by a transformation from the equivalence group SL(5). 

(b) System belongs to Monge-Ampere class (23). 

(c) System is invariant under an 8-dimensional subgroup of SL(5). 

(d) The characteristic variety defines conformal structure that is flat on every solution. 

Proof: 

equivalence (a) ( b ): One needs to show that taking a non-degenerate linear system, say 

u 3 = ^i, v 3 = u 2 (all non-degenerate linear systems of type (1) are SL(5)-equivalent), and 
applying transformations from the equivalence group, one obtains all systems of Monge-Ampere 
type. The easiest way to see this is the following. First of all, any non-degenerate linear system 
is invariant under an 8 -dimensional subgroup of SL(5); note that 8 is the maximal possible value 
for the dimension of the stabiliser of the action of the equivalence group on the space of fourfolds 
X C Gr(3,5). The corresponding infinitesimal generators can be calculated using the standard 
machinery of group analysis. Thus, the linear system u 3 = v 1 , v 3 = u 2 is invariant under the 
subgroup with 8 infinitesimal generators, 

U l5 V 2 , U 2 + V 3 , Ug + Vi, 

X13 + 2 X 32 + L12, X 23 + 2 X 3 i + L 2 i, Xn+X 22 + X 33 , Xu — X 22 + Lu, 

(we use the notation of Sect. 3). The Lie algebra is isomorphic to a semi-direct product 
V 3 x 0 l( 2 ) where V 3 ~ 5 3 (M 2 ) is the irreducible 4-dimensional representation of sl(2) (and 
hence of fll(2) = sl(2) © R, with non-trivial action of R). Thus, applying to a linear system 
transformations from the equivalence group, one obtains a variety of systems depending on 
24 — 8 = 16 essential parameters. It remains to point out that the class of Monge-Ampere 
systems also depends on 16 essential parameters. 

equivalence (a) (c): The first implication follows from the fact that any non-degenerate 

linear system is invariant under an 8 -dimensional subalgebra of SL(5). To establish the converse, 
let G be the symmetry group of system ( 1 ). We can always assume that the point o, specified 
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by Ui = Vi = 0, belongs to the fourfold X C Gr(3, 5) corresponding to our system. Let G 0 be 
the stabiliser of this point in G. Note that dim G — dim G a < 4, as G takes X to itself. The 
stabiliser P of the point o is spanned by infinitesimal generators X,j, L^, Pj, Q.;. Since the 
system is non-degenerate, we can always bring it to a canonical form: 

u 3 = vi + o(ui, Vi), v 3 = u 2 + o{ui, Vi). (24) 

This form (together with the point o) is stabilised by the following four elements of P: 

X13 + 2X32 + L12, X23 + 2X31 + L21, X11+X22 + X33, Xn-Xaa + Ln. 

Thus, dim G 0 < 4 so that dimG < 8 . The equality holds only if dimG 0 = 4. However, the 
generator Xu + X 22 + X 33 acts by non-trivial rescalings on terms of order 2 and higher in (24). 
Hence, for dim G a = 4, all higher-order terms must vanish identically, leading to a linear system. 

equivalence (a) (d): Representing system (1) in evolutionary form ( 8 ) and introducing 

the corresponding conformal structure g (see the proof of Theorem 1 for explicit formulae), the 
condition responsible for conformal flatness in three dimensions is the vanishing of the Cotton 
tensor, 

\7 r (Rpq -^Rgpq) — Vq(R pr jRg pr ), (25) 

where R pq is the Ricci tensor, R is the scalar curvature, and V denotes covariant differentiation in 
the Levi-Civita connection of g. Calculating (25) and using ( 8 ) and its differential consequences 
to eliminate all higher-order partial derivatives of u and v containing differentiation by x 3 , we 
obtain expressions that have to vanish identically in the remaining higher-order derivatives (no 
more than fourth-order derivatives will occur in this calculation). In particular, equating to zero 
coefficients at fourth-order derivatives of u and v, we obtain the system of second-order PDEs 
for / and g: 

t _ 2 9u, r 

JuiUi g Vi —fu i ^ UiVi ' > 

f _ 9u i t 1 9u i f 

J u i u j g Vi -fuG UiVi ' gvi-fu, Ju i v i' 

1 b J J 

r I r _ f u j ~ 9v j j- _|_ fuj ~9vj j? 

JuiVj ~r JujVi — f Ui —g Vi JuiVi ~r f u .-g v . JUjVji 

here i,j take any values from 1 to 2 ; the equations for g can be obtained by the simultaneous 
substitution / «->• g and u ■<->• v. The general solution of system (26) leads to Monge-Ampere 
systems. Furthermore, modulo (26) all remaining components of the Cotton tensor vanish 
identically (in fact, this follows from the linearisability of Monge-Ampere systems in 3D). This 
finishes the proof of Proposition 3. 

Remark. Equations of Monge-Ampere type have a clear geometric interpretation. Recall that 
the Grassmannian Gr(3, 5) is (locally) identified with the space of 2 x 3 matrices U. Minors 
of U define the Pliicker embedding of Gr(3, 5) into projective space IP 9 . We identify Gr(3, 5) 
with the image of this embedding, which is a non-singular algebraic variety of degree five. Thus, 
Monge-Ampere systems correspond to sections of Gr(3, 5) by subspaces P' C P 9 . 

We emphasize that the linearisability of Monge-Ampere systems is an essentially three- 
dimensional phenomenon: in higher dimensions there exist integrable non-linearisable examples 
of Monge-Ampere type, e.g. 

V >2 — V \ = 0, U3U4 — U4U3 — 1 = 0, 


Vv 


f — _ L _ f 

J y iVi ~ f u .-g v . JUiVi 


fviUj fu i -gv i f u i v i "b f u —g v . f u i v j 5 (26) 
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which is equivalent to the Plebanski first heavenly equation [49]. 


4.2 Linearly degenerate systems 

From the point of view of their ‘complexity’, linearly degenerate systems come next after lin- 
earisable systems. In this section we obtain a complete list of normal forms of linearly degenerate 
integrable systems: all of them come from quadratic maps P 4 ---> Gr(3,5) given by the classical 
construction of Chasles. 

The definition of linear degeneracy is inductive: a multi-dimensional system is said to be 
linearly degenerate (completely exceptional [6]) if all its traveling wave reductions to two dimen¬ 
sions are linearly degenerate. Thus, it is sufficient to define this concept in 2D case, 


U 2 = v 2 = g(ui,vi). (27) 

Setting u\ = a, v\ = p and differentiating by x l one can rewrite this system in two-component 
quasilinear form, 

«2 = /(o,p)i, P 2 = g{a,p)i, 

or, in matrix notation, 

f a \ . f a \ . f f a f p \ 

UA^UA- A= (l %)■ 

Recall that the matrix A is said to be linearly degenerate if its eigenvalues (assumed real and 
distinct) are constant in the direction of the corresponding eigenvectors. Explicitly, L r i A* = 0, 
no summation, where L r i denotes Lie derivative in the direction of the eigenvector r l , and 
Ar l = AV*. For quasilinear systems, the property of linear degeneracy is known to be related, 
under appropriate ‘small norm’ assumptions, to the impossibility of breakdown of smooth initial 
data, leading to global solvability of a generic Cauchy problem [52]. In terms of the original 
functions f{u\,v{) and g(u\,vi), the conditions of linear degeneracy reduce to a pair of second- 
order differential constraints, 


(/«! 


gvi)fuiui T 2 gmfuivi T gu\Pv\v\ T fv^guiui — 0, 

f ui )g V\V\ T 2 fviPui vi + t guifvivi — o. 


(28) 


Equations (27) define a 2-parameter family of 2-dimensional linear subspaces of a vector space 
V 4 (on projectivisation, a congruence of lines in P 3 ). Geometrically, equations (28) mean that 
focal surfaces of this congruence degenerate into curves [2]. We point out that system (28) 
is invariant under the natural action of the 2D equivalence group SL(4). Requiring that all 
traveling wave reductions of a multi-dimensional system to 2D are linearly degenerate in the 
above sense, we obtain differential characterisation of linear degeneracy: 

Proposition 4. A d-dimensional system, Ud = f(ui,vi), vj = g(ui,Vi), i = l,...,d — 1, is 
linearly degenerate if and only if the functions f and g satisfy the relations 

S v m {i,j,k } ((./Afc — gv k )fuiUj + g Uk (fu iVj + fujVi ) + fv k gu iUj + gu k gviVj) = o, 

(29) 

S V m {i,j,k} ((9 v k fu k )gv iVj + fv k (guiVj + g U jVi ) + gu k fviVj + fv k fmuj) — o, 

where Sym denotes complete symmetrisation over i,j,k€{l,...,d — l}. 
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Proof: 


Let us begin with the 3-dimensional case, 

u 3 = f(ui,u 2 ,v 1 ,v 2 ), v 3 = g(u 1 ,u 2 ,v 1 ,v 2 ). 


Looking for a traveling wave reduction in the form 

u^x^x 3 ) = u(y\y 2 ) + y 3 , u(x\x 2 ,x 3 ) = u(y\y 2 ) + y 4 , 

where y l are arbitrary constant-coefficient linear forms in the independent variables x k , one 
obtains a system of type (27) for u(y 4 ,y 2 ) and v(y 1 ,y 2 ). It is required to be linearly degenerate 
for every choice of linear forms y l . This condition imposes a set of second-order differential 
constraints for / and g, 


( fuj 9vj) fuiUi T 2( fm 9vi)fuiUj T 2 gujfmvi + 2 g UiifuiVj T 
fvjduiUi ffi 2 fviduiUj T 9uj9viVi T 2 9m9viVj = 0, 

(9vj ~~ fuj)9viVi T 2 (g Vi fui)9viVj + 2 fvj9uiVi + 2 fvi{9uiVj T 9ujVi)~\~ 

9ujfviVi + 2 gUifviVj T fvj fuiUi T 2 fvifuiUj = 0, 


(30) 


(no summation), where i,j £ {1,2} (8 relations altogether). Notice that the second half of these 
relations follows from the first one under the simultaneous substitution / -H- y, u •£)• v. Relations 
(30) can be obtained from (29) by setting k = j. For i = j = 1 relations (30) simplify to (28). 
Relations (30) can be derived in either of the two ways: 

Method 1. Applying the full 3D equivalence group SL(5) to constraints (28) one gets 8 linearly 
independent relations that are equivalent to (30). One can show that the set of relations (30) 
is SL(5)-invariant. This procedure generalises to higher dimensions in an obvious way. For 
instance, in 4D, applying transformations from SL(6) to constraints (28) one gets 20 linearly 
independent relations (29), d = 4. 

Method 2. Looking for particular traveling wave reductions of the form ^(x 1 ,^ 2 ,® 3 ) = 
u {y 1 -,y 2 )i u(x\x 2 ,x 3 ) = v(y 1 ,y 2 ) where y 1 = x 1 + Ax 2 , y 2 = x 3 , one obtains the reduced 
system in the form 


Uy2 F^Uyl^Vyl) f ^Uyl , Xllyl , Vyl , \V yl ) , G(Uyl , Vy\ ) g (lL yl ^ \U yl ^ V y\ ^ \V yl^j . 

Here F and G have to satisfy relations (28). Since F u t = f Ul + F Uy lU y l — fuiUi + 2A f Ul u 2 + 
fu 2 u 2 i etc, both relations (28) will be polynomial of degree 3 in A. Equating to zero coefficients 
of these polynomials we get all relations (30). This generalises to higher dimensions: thus, in 
4D, given a system 


ua = f(ui,u 2 ,u 3 ,v 1 ,v 2 ,v 3 ), va = g{ui, u 2l u 3 ,v ll v 2l v 3 ), 

one looks for traveling wave reductions ^(x 1 ,x 2 ,x 3 ,x 4 ) = ^(y 1 ,?/ 2 ), u(x 4 ,x 2 ,x 3 ,x 4 ) = v(y 1 ,y 2 ) 
where y 1 = x 1 + Ax 2 + /rx 3 , y 2 = x 4 . For the reduced system, both relations (28) become 
polynomials of degree 3 in A and /i, each having 10 coefficients. Equating them to zero we 
obtain 20 differential constraints (29) constituting conditions of linear degeneracy in 4D. This 
finishes the proof of Proposition 4. 
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Remark 1. Second-order relations (30) governing linearly degenerate systems in 3D are not 
in involution, and their prolongation implies all third-order integrability conditions (16). This 
requires differentiating equations (30) three times and solving for the fifth, fourth and third-order 
partial derivatives of / and g in terms of the first and second-order derivatives. 

Proposition 5. The parameter space of linearly degenerate integrable systems in 3D is 22- 
dimensional. Furthermore, the following conditions are equivalent: 

(a) System (1) is linearly degenerate and integrable. 

(b) There exists a unique flat symmetric connection on X in which the associated GL(2,M) 
structure is parallel. 

Proof: 

One can verify that 8 constraints (30) are compatible with integrability conditions (16). Re¬ 
membering that the parameter space of integrable systems (1) is 30-dimensional, we thus obtain 
a 30 — 8 = 22-dimensional parameter space of 3D linearly degenerate integrable systems. 

The equivalence of (a) and (b) can be demonstrated as follows. Using parametric equations 
of the fourfold X as in Theorem 1, 

u x = a , u y = b, v x =p, v y = q, u t = f(a,b,p,q), v t = g(a,b,p,q), 

we define cubic cones of the induced GL(2,R) structure as the intersection of three quadratic 
forms, u a = 0 , where 

oj 1 = dadq — dbdp, J 2 = dadg — dpdf , cj 3 = dbdg — dqdf. 

Note that these quadratic forms are nothing but restrictions to X of the second fundamental 
forms of the Grassmannian Gr(3,5) C P 9 . A connection V preserving this GL(2,R) structure 
is defined by the relations Vw“ = 0 mod(w^). Under the assumption that V is symmetric, 
these relations lead to a linear inhomogeneous system for the 40 Christoffel’s symbols of V. 
This system is uniquely solvable if and only if the functions / and g satisfy conditions of linear 
degeneracy (30). One can verify by direct calculation that the vanishing of the curvature of this 
symmetric connection is equivalent to integrability conditions (16). This finishes the proof of 
Proposition 5. 

Our next goal is to provide a complete list of normal forms of linearly degenerate integrable 
systems. The key example generating an open part of the 22-dimensional parameter space is as 
follows: 

Example. Let 01 , 02,03 and 01 , 02,03 be two triplets of numbers such that ai + 02 + 03 = 0 
and di + 02 + a 3 = 0 . One can show that the following system, 

a\a2U x Vy a2a\UyV x — 0, a\a^u x vi a^a\UfV x — 0, (31) 

is linearly degenerate and integrable. Eliminating v we obtain a second-order equation for u, 
a\u x u y t + a 2 U y u x t + a^utu X y = 0. Similarly, eliminating u we obtain the analogous equation 
for v, d\v x v y t + CL 2 V y v x t + 03 vtv xy = 0. This construction first appeared in [60] in the context 
of bi-Hamiltonian systems and Veronese webs in 3D, see also [21]. Rewriting system (31) in 
simplified form, 

u x v y = au y v x , u x v t = (3u t v x , 
a, /3 = const , one can show that: 
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(i) Systems with generic values of a and j3 are not SL(5)-equivalent. 

(ii) For generic choice of a and (3 the corresponding system is invariant under a 4-dimensional 
subgroup of the equivalence group (with infinitesimal generators Xu, X 22 , X 33 , Lu in the 
notation of Sect. 3, that form a Cartan subalgebra of SL(5)), and hence generates a 24 —4 = 20- 
dimensional orbit. 

Acting on (a, /3)-systems by transformations from the equivalence group one thus obtains a 
22 -parameter family generating an open part of the parameter space of linearly degenerate 
integrable systems. 

All linearly degenerate integrable systems in 3D, including the above example, can be ob¬ 
tained from the following geometric construction. Consider projective space P 4 with homoge¬ 
neous coordinates £ = (£ 4 : £ 2 : £ 3 : £ 4 : £ 5 ). Let A be a projective automorphism of P 4 defined 
by a 5 x 5 matrix from SL(5). Consider the family of lines (£, 77) through £ and 77 = £A (the locus 
of lines spanned by an argument and the value of a projective transformation was apparently 
first discussed by Chasles [13]; see also [14], p. 556). The Pliicker coordinates p lJ = 
are the 2 x2 minors of the matrix 


/ £2 £5 

^ 77 1 77 2 77 s 77 4 77 s 


they define quadratic map from the family of lines (£, 77 ) (which is itself isomorphic to P 4 ) into 
the Grassmannian of lines in P , that is, into Gr(2,5). By duality, this gives a fourfold X 
in Gr(3, 5), and the corresponding system (1). Fourfolds X arising from this construction are 
images of quadratic maps P 4 > Gr(3, 5) C P 9 . Explicit parametric equations of X , as well as 
of the corresponding system ( 1 ), can be obtained from the factorised representation 


e 1 £ 2 £ 3 £' £ 5 

1 0 S 4 K 

77 77 77 77 77 


£' £ 2 

1 2 

77 77 


explicitly, 


Ul = 


P 


32 


pl2 p±* pit, p 

In the (generic) diagonal case, A = diag (\ l ), this gives 


u 2 = 


P 


,42 


,12 


, U 3 = 


P' 


52 


,12 : 


Vl = 


P 


1 0 «2 7X3 

0 1 v\ v 2 v 3 


13 


12 ■ 


V2 = 


14 


12 ' 


pl5 

v 3 = 12 ' 

pLZ 


A 2 - A 3 £ 2 £ 3 

A 2 - A 4 £ 2 £ 4 

A 2 - A 5 £ 2 £ 5 

Ul ~ A 2 — A 1 £ 4 £ 2 ’ 

U2 ~ A 2 - A 1 £ 4 £ 2 ’ 

U3 ~ A 2 - A 1 £ 4 £ 2 

A 3 - A 1 £ 4 £ 3 

A 4 - A 1 £ 4 £ 4 

A 5 - A 1 £ 4 £ 5 

Vl ~ A 2 -A 1 £ 4 £ 2 ’ 

V2 ~ A 2 - A 1 £ 4 £ 2 ’ 

“ A 2 -A 1 £ 1 £ 2 ' 


leading to the (a,/3)-system, u\v 2 = au 2 v±, u±v 3 = (3u 3 v\, with 


(32) 


(A 2 -A 4 )(A 3 -A 4 ) (A 2 -A 5 )(A 3 -A 4 ) 

(A 2 - A 3 )(A 4 - A 1 ) ’ ° “ (A 2 — A 3 )(A 5 — A 1 ) ’ 


The corresponding fourfold X coincides with the quadratic image of P 4 defined by a linear system 
of quadrics that pass through 5 points in general position (in our parametrisation, quadrics 
£*£ J , i / j, that pass through 5 base points each having only one nonzero coordinate £). We 
point out that smooth quadratic embeddings of projective spaces into Grassmannians were 
discussed in [54], It is important to emphasize that our X is not smooth: lines through any 
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pair of base points in P 4 correspond to singular points of X. In algebro-geometric language, 
quadratic maps from P" to the Grassmannian of lines in P n are given by rank 2 vector bundles 
E with the first Chern class equal to 2, together with n+1 sections that generate E. In the 
Chasles construction, E = 0(1) + 0( 1), and the space of sections is the graph of the map 
A : H 0 (O( 1)) -> H 0 (O( 1)). 

Table 1 below comprises a complete list of 7 canonical forms of linearly degenerate integrable 
systems (1) labelled by Segre types (essentially, Jordan normal forms) of 5 x 5 matrices A from 
the Chasles construction. Thus, type [3, 2] denotes two Jordan blocks of size 3x3 and 2x2, 
type [5] - one 5x5 Jordan block, etc (to get canonical forms as presented in Table 1 we take A 
in Jordan normal form, eliminate £’s from (32), and use equivalence transformations to remove 
unessential parameters). The last column provides dimensions of stabilisers of the resulting 
systems under the action of the equivalence group SL(5). 


Table 1: Canonical forms of linearly degenerate integrable systems in 3D 


Segre type 

Canonical form 

Dim(stab) 

[11111] 

UlV2 = OLU2VI 

U1V3 = PU3V1 

4 

[2111] 

U1V2 — U2V1 = V1V2 

U1V3 — U3V1 = av 1V3 

4 

[221] 

U1V2 — U2V1 = v\ 

U1V3 — U3V1 = V1V3 

4 

[311] 

IL2 = V1V2 

U3 = (1 — V\)V3 

4 

[32] 

U2 = V1V2 

U3 = V2 + V1V3 

5 

[41] 

Ul = V2 — Vi 

U3 = (1 — Vl)v3 

5 

[5] 

Ul = V2 — v\ 

U2 = V3 — V\V2 

6 


Systems presented in Table 1 are not SL(5)-equivalent. Even though all 4-dimensional stabilisers 
are commutative subalgebras of SL(5), they are not conjugate (in particular, only the first of 
them is a Cartan subalgebra). 

Remark 2. Note that Table 1 does not contain a linear system. It can be recovered from the 
following generalisation of the Chasles construction. Take projective space P 4 with homogeneous 
coordinates £ = (^ 4 : £ 2 : £ 3 : £ 4 : £ 5 ), let A and B be two 5x5 matrices from SL(5). Consider 
the family of lines (rj, (p) through 7] = £A and = £B. A linear system corresponds to the case 
when the pencil A + A B consist of two Kronecker blocks of the size 1x2 and 4x3 (Kronecker 
blocks of other size lead to degenerate systems). 

Proposition 6 . The Chasles construction gives all linearly degenerate integrable systems in 

3D. 


Proof: 


Consider the Grassmann variety of all pencils A+XB in the space of 5 x 5 matrices. Denote by 
£ the open subset in Gr(2, 25) which corresponds to non-degenerate equations via the Chasles 


23 













construction. Clearly, this is an irreducible algebraic variety. The Chasles construction defines 
a rational algebraic map from Gr(2,25) to the 22-dimensional variety of all linearly degenerate 
integrable systems, with £ as the set of regular points. The fibres of this map are orbits of the 
right action of SL(5) on pencils A + A B. Moreover, we can easily verify that this action has a 
discrete stabiliser in the case when the pencil A + \B is of Segre type [11111]. Therefore, generic 
fibres of this map are 24-dimensional. Thus, its image is an irreducible variety of dimension 
2 ■ 23 — 24 = 22. As Gr(2,25) is a projective variety, the image has to coincide with the set of 
all linearly degenerate integrable systems. 

Table 2 brings together some further examples of linearly degenerate integrable systems that 
can be interpreted as Backlund transformations: on elimination of v, they lead to a second-order 
PDE for u, similarly, on elimination of u they lead to a second-order PDE for v (columns 2 and 
3). Each of these systems is equivalent to one of the canonical forms presented in Table 1 (as 
indicated by its type). 


Table 2: Examples of Backlund transformations in 3D 


System (1) 

Equation for u 

Equation for v 

Type [11111] 

aia,3Vtu x — a,3a\v x ut = 0 

aia 2 VyU x — a, 2 aiv x u y = 0 

CL\UxUyt + U^UyUxt U^UtUxy — 0 

fii v x v y t + 0 , 2 VyVxt + a,3VtV xy = 0 

Type [2111] 

(A UyVx - 0 

A vt — u t v x = 0 

Uyt + UyUxt UtUxy — 0 

VxVyt + (A - 1 )VyVxt — A VtVxy = 0 

Type [221] 

Uy \UyVx - 0 

Vt + {^ 2 u y — \u t )v x = 0 

Uyy "b UyUxt Ut'U'X'y = 0 

UtUxy UyUxt “b X{UyUxy UxUyy ) = 0 

Type [221] 

UyVt — 1 = 0 

U X UyVx = 0 

'U'Xy “b UyUxt UxUyt — 0 

Vxt + VtVxy ~ VxVyt = 0 

Type [311] 

Vt + (A — Uy)Vy = 0 
\v x U X Vy = 0 

Uxt “b UxUyy UyUxy = 0 

UtUxy UyUxt ~b X{UyUxy UxUyy ) = 0 

Type [32] 

UyVx — 1=0 

Ut ~ UyVy = 0 

Uyy + UyUxt UtUxy — 0 

Uxt “b UxUyy UyUxy = 0 

Type [41] 

Vy + (X — u x )v x = 0 

vt + (A 2 — \u x + u y )v x = 0 

Uxt “b Uyy "b UyUxX UxUxy = 0 

(::)„ = 0v), 


Some of the second-order PDEs from Table 2 have appeared in different contexts in [47, 42, 60, 
10, 45, 3]. 

Remark 3. The equivalence group acts algebraically on the parameter space of linearly generate 
integrable systems, so we can consider the rational quotient M r 2 C g/SL(5). Since generic orbits 
of this action have dimension 20, the moduli space of linearly degenerate integrable systems is 
two-dimensional. 


24 













4.3 Further examples in terms of modular forms and theta functions 

This section contains a list of more exotic integrable systems that are not expressible in elemen¬ 
tary functions. These examples demonstrate that one should not expect any ‘simple parametri- 
sation’ of integrable systems of type ( 1 ). 

Example 1. Let us begin with the system 

Uy 1 . \ 2 

Ut = V X , Vt = — + -v(u x )v x , 
v x 6 

which appeared in [46] in the classification of integrable hydrodynamic chains. In this case the 
integrability conditions result in the Chazy equation for rj, 

rf" + 2rjr," - 3r/ 2 = 0, 


whose generic solution is known to be the weight two Eisenstein series E 2 associated with the 
full modular group. 

Example 2. Consider the system 

v x + u x u y r(ut ) = 0 , v t = u y . 

In this case the integrability conditions result in the third-order ODE for r, 

r'"{r' — r 2 ) — r" 2 + 4 r 3 r" + 2r /3 — 6r 2 r /2 = 0, 


which appeared recently in different context in the theory of modular forms of level two: compare 
with equation (4.7) from [1], Its generic solution is given by the Eisenstein series 


r{u t ) 


i-*E 

n=1 


(-1 ) n ne 4nut 
\ _ g4 nut ’ 


associated with the congruence subgroup To( 2 ) of the modular group. 
Example 3. As a generalisation of Example 2, let us consider the system 


v x + u x r(u y , u t ) = 0 , v t = u y . 

In a somewhat different representation, it first appeared in [10]. It was demonstrated that 
the integrability conditions imply r = 2 s'/s where prime denotes differentiation by u y , and 
s(u y ,ut ) = 0 ( 7 ^, — *4) is the Jacobi theta function: 

OO 

9(z, t) = 1 + 2 e mn2l ~ cos(27 rnz). 

n —1 


Example 4. Further generalisation, 


v x + f(u x ,u y ,ut) = 0 , v t = u y , 


was discussed in [25], where it was shown that the requirement of integrability implies that 
generic / is given by the ratio of two Jacobi theta functions: 


f(u x ,u y ,ut ) 


1 6i(u t ,u y - u x ) 
4 0 1 (u t ,uy + u x )' 


25 




4.4 Integrable systems in terms of generalised hypergeometric functions 

In this section we utilise the construction of [44, 45] which parametrises integrable quasilinear 
equations in 2 + 1 dimensions in terms of solutions of the generalised hyper geometric system, 

d 2 h Si dh Sj dh . 

a a = a ^ "a ’ 1 ' ■?> 

OZiOZj Zi — Zj OZj Zj — Zi OZi 


d 2 h 

8z 2 



72+2 

+E‘ 

3 = 1 


Zi(Zi - 1 ) 


h + 


Sj 

Zi(Zi - 1 ) 


n 


E 


*j(*j ~ x ) 

Zj - Zi 


dh 

dzj + 


n 



s i + 
Zi 


Si + S n +2 
Zi - 1 


dh 

<9+ 


(33) 


here h is a function of n variables zi,...,z n , and si,...,s n +2 are arbitrary constants. This 
system is involutive and possesses n + 1 linearly independent solutions known as generalised 
hypergeometric functions [29, 44], Here we consider the case n = 4. Let {h 1 , h 2 , <+ g 2 , g 3 } be a 
basis of solutions of system (33). Let us introduce the parametric formulae 


W{g\h}) 

W{g 2 ,h l ) 

W(g 3 ,h l ) 

Ux IT(h 1 ,/i 2 )’ 

dy W{h\h 2 )’ 

Ut w(h i ,h 2 y 

(34) 

W(g\h 2 ) 

W(g 2 ,h 2 ) 

W(g 3 ,h 2 ) 

Vx W(/i 1 , h 2 ) ’ 

Vy W{h l , h 2 ) ’ 

Vt W(h 1 ,/i 2 )’ 

where W is the Wronskian, W (/, g) = 

fg' — gf , and prime denotes differentiation with respect 


to one of the variables z* (for definiteness, we assume ' = d/dz 4). Formulae (34) can be viewed 
as parametric equations of a fourfold X C Gr(3, 5). Eliminating the parameters z\...., Z4 we 
obtain system of type (1). 

Proposition 7. Parametric equations (34) define generic integrable system of type (1). 

Proof: 

Given hypergeometric system (33), the construction of [44] requires an additional ingredient, 
namely, the choice of a fc-dimensional subspace of solutions of system (33). We will only need 
a particular case of the general scheme that corresponds to n = 4, k = 2. Thus, let h 1 ,/? 2 
be a pair of linearly independent solutions of system (33). Let us complete this pair to a 
basis {h 1 , h 2 , g 1 , g 2 , g 3 }, and introduce the following system of 6 first-order quasilinear PDEs for 
Zi ,..., Z4, considered as functions of the auxiliary variables x, y, t: 


W(g 1 ,h 1 ) \ _ ( W(g 2 y) \ 

W(hl,h*)) y - \W (h 1 ,h' 2 ) J x ’ 

W(g 1 ,h 2 )\ _ fW{g 2 ,h 2 )\ 

WW7?j) y - \ W(h L ,h 2 )) x ’ 


f W(g 1 ,h 1 ) \ _ ( W(g 3 ,h 1 ) \ 

\ W (h L ,h 2 ) J t ~ \W(h\h 2 )) x > 

fW(g 1 ,h 2 )\ _ fW(g 3 ,h 2 )\ 

\W(h\h 2 )) t - \ W(h\h 2 )) x > 


fW(g 2 ,h 1 )\ _ fW(g 3 ,h 1 )\ 

\W(h\h 2 )J t ~ \W(h\h 2 )J y » 

(W{g 2 ,h 2 )\ _ fW(g 3 ,h 2 )\ 

\W(h\h 2 )J t - y W (h\K 2 )) y - 


It was demonstrated in [44] that this system is integrable by the method of hydrodynamic reduc¬ 
tions, and possesses a dispersionless Lax pair. The conservative structure of this system implies 
the existence of potentials u and v specified by (34). Thus, parametric equations (34) indeed give 
rise to integrable system (1). Since the moduli space of integrable systems (1) is 6-dimensional 
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and, for n = 4, hypergeometric system (33) depends on 6 essential parameters si,... ,sq, the 
claim follows. One can show that transformations from the equivalence group SL(5) are in 
one-to-one correspondence with linear transformations of the basis {h 1 , h?, g 1 , g 2 , g 3 }. 

We emphasise that, although formulae (34) parametrise generic integrable systems (1), the 
degeneration procedure leading to examples from Section 4.3 is far from trivial. 

5 Geometric aspects of integrability 

In this section we adopt a geometric point of view and consider 3D system (1) as equations 
defining a fourfold X in the Grassmannian Gr(3, 5). Our aim is to reformulate the integrability 
conditions in intrinsic geometric terms. 

Recall that in Sect. 3 we introduced the main object needed to develop differential geometry of 
submanifolds of Gr(3, 5), namely, the field of Segre cones in the tangent bundle TGr(3, 5). The 
field of Serge cones is invariant under the natural action of the equivalence group SL(5), and 
plays the role of an ambient flat conformal structure. 

In Sect. 5.1 we develop geometry of fourfolds X C Gr(3, 5). The induced field of Segre cones 
gives rise to a field of twisted cubics in the projectivised tangent bundle PTX. Thus, a fourfold 
in Gr(3,5) carries an intrinsic GL(2,R) structure. 

In Sect. 5.2 we show that two-component and three-component hydrodynamic reductions of 
system (1) correspond to bisecant and holonomic trisecant submanifolds of the fourfold X, 
respectively. This allows us to reformulate the integrability geometrically as the existence of an 
infinity of holonomic trisecant submanifolds parametrised by three arbitrary functions of one 
variable. 

In Sect. 5.3 we demonstrate that every fourfold X corresponding to integrable system (1) carries 
a canonical conformally symplectic structure. 

In Sect. 5.4 we reformulate the integrability conditions in terms of the torsion/curvature of the 
Bryant connection of the associated GL(2,R) structure. 

5.1 GL(2,M) structure on a fourfold X C Gr(3,5) 

Let X be a non-degenerate fourfold in the Grassmannian Gr(3, 5). Taking a point o £ X and 
projectivising the intersection of the tangent space T 0 X with the Serge cone C in T 0 Gr(3,5) one 
obtains a twisted cubic, that is, a rational normal curve of degree three. Indeed, the projectivised 
Serge cone is a non-singular rational cubic threefold in P 5 , so that its section by the projectivised 
tangent space PT 0 A = P 3 is a non-singular curve of degree three, that is, a twisted cubic. This 
twisted cubic can also be interpreted as the set of matrices of rank one in the tangent space 
T 0 X, when we identify T 0 Gr(3, 5) with the space of 2 x 3 matrices. Thus, the projectivised 
tangent bundle of the fourfold X is equipped with a field of twisted cubics. This supplies X 
with a GL(2,R) structure. 

Alternatively, one can say that each tangent space to X is identified with a four-dimensional 
space of binary cubics: in this picture the rational normal curve corresponds to cubics with a 
triple root. In representation-theoretic language, the tangent space T a X is an irreducible sl(2) 
module V 3 ; here and below Vi — ^(R 2 ) is the irreducible sI( 2 )-representation of dimension l + 1 . 

Remark 1. Another natural source of GL(2, R) structures in four dimensions is provided by 4- 
th order ODEs with vanishing Wiinschmann invariants [ 8 , 15, 20, 43]. In this context, GL(2,R) 


27 


structures are induced on the parameter spaces of solutions to the corresponding ODEs. We 
emphasise that the structures coming from ODEs always have zero torsion (so that all geometry 
is contained in the curvature), while in our case the torsion is generally nonzero, and defines the 
curvature. 

Remark 2. An interesting class of GL(2, M) structures in five dimensions appeared recently 
in the context of integrable equations of the dispersionless Hirota type. These structures are 
induced on hypersurfaces of the Lagrangian Grassmannian that correspond to such equations, 
see [26, 55] for further details. 

The most important property of GL(2,R) structures in four dimensions is the existence 
of a unique intrinsic connection (the Bryant connection) whose torsion takes values in the 8 - 
dimensional irreducible representation of GL(2, R) [ 8 ]. Let us recall the construction. The cubic 
cones of the GL(2,R) structure induced on X are defined as the intersection of 3 quadrics, 
ui a = 0 , where 

uj 1 = dadq — dbdp, oj 2 = dadg — dpdf, u; 3 = dbdg — dqdf ; 

here f(a,b,p,q ) and g(a,b,p,q ) correspond to the right-hand sides of system ( 8 ). The tangent 
space T 0 X is an irreducible sl(2) module V 3 , and since A 2 V 3 = Vo + V 4 , it has a canonical 
conformal almost symplectic structure D (see Sect. 5.3 for explicit formulae). Raising one index 
of w“ by means of D we obtain 3 operators A a £ g 1(4), a = 1,2,3, that form a basis of the 
irreducible embedding of sl(2) into g[(4). Define the bilinear form 

B ap = tiA a A^ 

which, up to a constant factor, is the Killing form of sl(2). Let B # be the inverse of B, and let 

C = 20B*pA a AP 

be the (normalized) Casimir operator. Then C is a well-defined element of the universal envelop¬ 
ing algebra of sl( 2 ) which is independent of the choice of the basis cj" (the coefficient 20 appears 
here because we use the embedding of sl(2) into flt(4), rather than the adjoint representation, 
to define the Casimir operator). The Casimir C acts naturally on any irreducible sl(2) module 
Vi as a scalar operator (due to the Schur lemma), 

C Vl =Xild, Xi = 1(1 + 2). 

The Bryant connection V is uniquely defined by the following two properties [ 8 ]: 

• Vw“ = 0 mod (ujP). This means that V preserves the GL(2,R) structure. 

• The torsion T of the Bryant connection V lies in the unique 8 -dimensional submodule Vj 
of the sl(2)-representation in the space of all algebraic torsions. This can be expressed via 
the following linear condition: 

C -T = X 7 T = 63T. 

In more detail, the first condition gives 56 equations for the 64 Christoffel’s symbols T l - k of 
the connection V. Of these 48 are linearly independent, so that V exists with the freedom of 
64 — 48 = 16 arbitrary functions. To fix these functions we note that since the tangent space 
r = T a X is an irreducible sl(2) module V 3 , the torsion T of V decomposes as follows: 

AV <S> t = (Vq © Vi) ® V 3 = Vf ® 2V 3 © V 5 © V 7 . 
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We can change the connection by a g[(2)-equivariant gauge from Hom(r, gl(2)), which decom¬ 
poses as V 3 <S> (Vo + V 2 ) = V\ © 2 V 3 © V 5 . This removes the freedom of 16 arbitrary func¬ 
tions and yields a unique connection with torsion T 6 Vj. Note that the action of the op¬ 
erators A a on tensors is standard, for instance, for the torsion tensor one has ( A a ■ T)V = 
(A a )^T^ — (A a )fT^ — (A Q )°T*. This defines the action of the universal enveloping algebra, in 
particular, of the Casimir operator C. 

5.2 Geometric interpretation of the integrability conditions 

As explained in Sect. 5.1, each projectivised tangent space T 0 X carries a rational normal curve 
7 of degree three. 

Definition. 

(a) Bisecant surface is a two-dimensional submanifold X 2 C X whose projectivised tangent 
planes are bisecant lines of 7 . 

(b) Given a three-dimensional submanifold X 3 C X, each tangent space T 0 X 3 carries three 
distinguished directions, namely those corresponding to the three points of intersection of PT 0 X 3 
with 7 . These directions define a net on X 3 , we will call it the characteristic net. Holonomic 
submanifolds X 3 are defined by the requirement that the characteristic net is holonomic (that 
is, locally a coordinate net). 

Proposition 8. Bisecant surfaces and holonomic submanifolds X 3 correspond to two- and three- 
component hydrodynamic reductions of system (1). Furthermore, 

(i) every fourfold X possesses infinitely many bisecant surfaces parametrised by two arbitrary 
functions of one variable; 

(ii) a fourfold X corresponds to an integrable system if and only if it possesses infinitely many 
holonomic submanifolds X 3 parametrised by three arbitrary functions of one variable. Thus, the 
existence of holonomic submanifolds X 3 is a geometric characterisation of integrability. 

Proof: 

We follow the notation of Sect. 2.1. Let us represent our system in quasilinear form (9), and 
take an N-conrponent reduction specified by 

a = a{R 1 ,..., R N ), b = b(R\...,R N ), p = ^(f ? 1 ,..., R N ), q = q(R x , ..., R N ). 

The geometric image in the Grassmannian Gr(3, 5) is a submanifold X^ C X (it would be 
sufficient for our purposes to restrict to N = 2, 3), represented by the symmetric matrix 

u=( a b f ), (35) 

V P Q 9 J 

parametrised by R 1 , ...,R N . Using equations (10), (11) one obtains the following expression for 
the derivative of U with respect to R 1 , 


( dia pddia X l dia \ 

V dip pddip Xdip ) ' 


(36) 


Thus, diU is a matrix of rank one, so that the projectivisation of diU belongs to the rational 
normal curve 7 , and coordinates R l provide the characteristic net on X jV . For N = 2 we 
have a two-dimensional surface X 2 C X parametrised by f? 1 ,^ 2 . Since both d\U and S 2 U 
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have rank one, the surface X 2 is bisecant. As explained in Sect. 2.1, every system of type (1) 
(not necessarily integrable) possesses infinitely many two-component reductions parametrised by 
two arbitrary functions of one variable. Thus, every fourfold X of the Grassmannian Gr(3,5) 
possesses infinitely many bisecant surfaces. This establishes part (i) of the Proposition. 

Similarly, every three-component reduction corresponds to a holonomic submanifold X 3 C 
X parametrised by R ] . R 2 , R 3 . Since an integrable system possesses (by definition) infinitely 
many three-component reductions parametrised by three arbitrary functions of one variable, the 
corresponding fourfold X possesses infinitely many holonomic submanifolds X 3 . This establishes 
the first part of (ii). 

To finish the proof one needs to show that, conversely, bisecant surfaces (holonomic subman¬ 
ifolds X 3 ) of A correspond to two-component (three-component) reductions of the associated 
system. This can be demonstrated as follows. Let X 2 be a bisecant surface represented in form 
(35), referred to its characteristic net R , 1 , R ?. Thus, a, b,p, q,f, g are functions of R l . R 2 such 
that the rank of diU equals one, so that one can introduce parametrisation (36). Compatibility 
conditions for the equations dib = gddia and diq = pddip imply 


didja 


. d,a + d 3 a and d l d r p 

/P — p , 1 p l — pP 


d jX 


jdiP + 


djp 3 
P 


1°^ 


respectively. Similarly, compatibility conditions for the equations dif = 
imply 


didna = . . 

3 X - A* 


d ’ x ' 9, n + 


A* - X 3 

Subtracting these equations from each other we get 


and didjp = 


djX 
X - A* 


dip + 


Xdid and dig 
diX 


A* - X 


d jP . 


Xdip 


Sijdid + Sjidjd = 0 , Sijdip + Sjidjp = 0 , 

where Sij = djpL l /(gi — //) — dj\ l /(X — A*). The case Sij = Sji = 0 gives commutativity 
conditions (7), so that we recover all equations governing two-component reductions. Thus, 
bisecant surfaces indeed correspond to two-conrponent reductions. If one of the coefficients 
or Sji is nonzero, then the functions a and p must be functionally dependent, so that dp = cda 
for some coefficient c. But then we must also have dq = cdb and dg = cdf, so that the whole 
tangent space of the surface X 2 consists of rank one matrices. This, however, is not possible for 
non-degenerate systems ( 1 ). 

For holonomic submanifolds X 3 , considerations are essentially the same. Since the existence 
of three-component reductions is necessary and sufficient for the integrability, this finishes the 
proof of Proposition 8 . 


5.3 Canonical conformal almost symplectic structure 

Let P 3 be a projective space with homogeneous coordinates (x° : x 1 : x 2 : x 3 ). A rational normal 
curve 7 of degree three has the standard form 7 = (1 : t : t 2 : t 3 ). One can verify that, up to 
a scalar factor, there exists a unique 2-form, oj = dx° A dx 3 — 3 dx l A dx 2 , for which all tangent 
planes to the cone over 7 are Lagrangian, that is, 07 ( 7 , 7 ') = 0. The matrix of co is 

/ 0 0 0 1 \ 

0 0-30 

0300' 

\ -1 0 0 0 / 
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Since each projectivised tangent space of X contains a rational normal curve, the ‘curved’ version 
of the above construction supplies X with a nondegenerate 2-form 17 defined up to a conformal 
factor, that is, a conformal almost symplectic structure. 

The differential of 17 can be represented in the form di 7 = p A 17 where p is the so-called 
Lee form of 17 (in four dimensions such representation exists, and is unique). Since the rescaling 
17 —> /!7 gives rise to the transformation p —»• p + df / f, the differential dip is independent of the 
conformal factor /. In particular, the condition dp = 0 is (locally) equivalent to the existence 
of a function / such that /17 is symplectic (such structures 17 are also known as conformally 
symplectic). 

Proposition 9. The canonical structure 17 on the fourfold X corresponding to integrable system 
(1) is locally conformally symplectic. If H l (X, M) = 0, and X is non-degenerate, then there 
exists a symplectic structure 17 on X defined canonically up to multiplication by a constant. It 
is parallel with respect to the Bryant connection. 

Proof: 

Following the notation of Sect. 2.1 we parametrise X C Gr(3, 5) by local coordinates a,b,p,q 
in the form 

X = ( ° b 

\ p q g(a, b,p, q) ) ' 

The field of rational normal curves 7 in the projectivised tangent space of X is specified by the 
equations rk(dX) = 1, that is 


da dq — dbdp = 0, da dg — df dp = 0, db dg — df dq = 0, 


or, explicitly, 7 = (m : n : mt : nt ) where m = g b + (g q - f b )t - f q t 2 , n = -g a + (f a ~ 9 P )t + f p t 2 , 
and t is a parameter. In matrix notation we have 7 = (1 : t : t 2 : t 3 ) ■ A, where 


A 


1 9b ~9a 0 0 \ 

9q fb fa 9p 9b ~9a 
— fq fp 9q — fb fa — 9p 

V 0 0 -fq fp J 


Thus, the matrix of the 2-form 17 is given by 17 = . It remains to calculate the 

corresponding 1-form p, and verify that the condition dp = 0 holds identically modulo integra- 
bility conditions (16) satisfied by / and g. We point out that defl7 = 9/(detA) 2 , furthermore, 
the condition detA = 0 is equivalent to the reducibility of dispersion relation (12). Thus, our 
assumption of non-degeneracy of the characteristic variety is equivalent to non-degeneracy of 
17. Finally, one can verify by direct calculation that the associated symplectic structure, which 
is unique up to a constant factor, is parallel in the Bryant connection: this follows from the 
fact that both the connection and the symplectic structure are defined in a canonical way. This 
finishes the proof of Proposition 9. 


We emphasize that the condition dp = 0 is necessary, but not generally sufficient for inte- 
grability. For instance, one can show that this condition is satisfied for every system of the form 
ut = f(u x ) + v y , u y = v x , where the function / is arbitrary, whereas the integrability implies 
f" = 0. On the other hand, for systems of the form ut = v y , vt = u y + g(u x ), the condition 
dp = 0 implies g'"g' — 2g" 2 = 0, which is precisely the integrability condition. We skip details 
of these calculations. 
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5.4 Integrability conditions via differential invariants 

Let V be the Bryant connection of a GL(2, Restructure on X. Its basic invariant is the torsion 
T G V 7 . Other tensorial invariants include the curvature R, and the covariant derivative of 
torsion, XT. Let us give here explicit coordinate formulae that show our conventions (for non- 
symmetric connections): 

X k di = T“ k d a , X k dx l = —T i ak dx a , 7g = (Xidj - Vjdi, dx k ) = T% - T^, 

Ruj = ([Vi, Vj]di, dx k ) = ditfj - djTi + - T° u T h av 

{X7 r r\k _ f) 'T'k 1 p& r r a -pa rpk pa rpk 

Notice that both R and VT (which are third-order differential invariants of the corresponding 
systems (8), as L')) are second-order quantities in partial derivatives of / and g ), belong to 
the space r © r* © A 2 r* which has the following sl(2) decomposition into irreducibles (with 
multiplicities): 

T <g) T* <8) AV = 2Vo © 414 0 514 © 414 © 2Vs © V w . (37) 

Consequently, we can decompose the basic third-order invariants as follows: 

R = R( 0 ) + R(2) + R (© + R( 6), vr = vT( 4 ) + VT( 6 ) + VT( 8 ) + VT( 10 ), 

where the subscript (A;) indicates the weight, or equivalently the (multiple) submodule V k to 
which the component belongs. Indeed, the spaces where these tensors take values decompose 
into sl(2)-irreducibles as follows: 

01(2) © A V = Vo © 214 © 214 © Ve, r* © V 7 = V 4 © V 6 © V 8 © V w . 

For a tensor K , the condition K G V/ is equivalent to C ■ K = 1(1 + 2)17. This gives a decompo¬ 
sition of i? and VT into eigenspaces of (7. 

By Proposition 5 of Sect. 4.2, system (1) is linearly degenerate and integrable if and only if 
the associated Bryant connection is trivial: T = R = 0 (this is also a consequence of Theorem 
3 below). It turns out that in the general, not linearly degenerate case, the integrability is 
characterised by the condition that both R and VT can be represented as certain quadratic 
expressions in the torsion T. To write down these expressions we introduce the following tensors 
T 2 , T 2 , T 2 , T 2 , T 2 G t©t*©A 2 t* (the upper index 2 indicates that these tensors are quadratic 
in T, and a, (3,^,5 are labels, not indices): 

(rp2 \ k rpk rpa ( rp2\k rpk (XdbrpC r\ ( rp2\k rpk (XClbrpC (") 

\-* )lij \-*-(3)lij ^\jcr^ 

(rji2 \ k rjik (XClbrjiC o frj-\2\k (xkarjib 0 rjic . 

— 1 [io} 1 1 bjY l cU (J-SJlij — -Lal^bc-Lij, 

here square brackets denote skew-symmetrization in i,j, and fl is the canonical almost symplectic 
structure (note that these tensors are independent of the choice of a conformal factor of II). Due 
to the identity T 2 = 2 T 2 , the above formulae give 4 essentially different invariant tensors that 
are quadratic in the second-order partial derivatives of / and g. For every a G {a, /3, 7 , 5} we 
decompose T 2 into sl(2)-irreducibles, where apriori all weights 0 < l < 10, l G 2Z, are possible: 

T '2 11 <2 I rp2 I rji2 | ' / '2 I rp2 1 rp2 

a ~ 1 a{ 0) + 1 a( 2) + 1 a ( 4) + 1 cr (6) + 1 ( r{ 8) + 

We claim that in fact 


r « 40 ) = 4) = T 2 (8) = 0 for all a G {a, P, 7,<5}- 
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Indeed, these tensors are contractions of T © T £ S 2 Vr = 14 © V§ © V 10 © V 14 , which is then 
projected to the right hand side of (37). Thus, Vj .4 disappears, and only components 14 with 
k = 2 , 6,10 remain in the decomposition. 

Theorem 3 For non-degenerate system (1), the integrability is equivalent to the following rela¬ 
tions among the invariants of the associated GL(2,R) structure: 

R(o) = 0) R(i) = 0, VT( 4 ) = 0 , VT (8) = 0 , VT( 10 ) = —28T^( 10 ), 

p _ 44 rp2 1 o nn2 40 n 

^(2) - T J a(2) + 2J /3(2) - — J 7 (2) ~ ^<5 (2), 

-^( 6 ) = — 24Tj( 6 ) - 30 T|( 6 ) — 60 ( 6 ) — 24 T|( 6 ), 

VT( 6 ) = - 8 Tj( 6 ) - 8 T|( 6 ) - 16( 6 ) - 4T|( 6 ). 

Note that for /c = 2,6 there are 4 linearly independent components of 

weight fc, and 14 has multiplicity 4 in the sl(2)-submodule = (r © r* © A 2 r*)( fc p Thus, 
it is expected that and VT4) decompose in this ‘basis’ (since the integrability conditions 
reduce third-order expressions to second-order). Similarly, for k = 10 just one T 2 ( 10 ) is a basis 
of n (10 ), and VT( 10 ) is expressed through it. However, the tensors T 2 have no 14-components 
for k = 0 ,4, 8 , and for the corresponding third-order tensors we get simpler relations. 

Proof: 

The proof is computational. To verify these formulae we substitute the integrability condi¬ 
tions (16). This leads to expressions that are quadratic in the second-order partial derivatives 
of / and g, with coefficients depending on their first-order derivatives (that is, functions on 
J 1 ). The computational complexity of the output is quite high, however, the substitution of 
any ‘generic’ point of J 1 readily gives zero (this verification is done in Mathematica; we choose 
a rational point so that all computations are exact and rigorous). Here the choice of a generic 
point is irrelevant due to the transitivity of SL(5)-action on J 1 , see discussion at the end of 
Sect. 3. 

Conversely, one has to verify that the above relations imply all of the 40 integrability condi¬ 
tions (16). Note that these relations are linear in third-order partial derivatives of / and g, with 
coefficients being functions on J 1 . The rank of the matrix at third-order derivatives is equal to 
40 (one can again restrict to any generic point of J 1 ), and this implies the claim. 

Remark. Some of the relations from Theorem 3 hold identically for every fourfold X C Gr(3, 5), 
without using the integrability conditions. These include the relations 

77) _ n _O nn2 rp2 _ rp2 rp2 _ rp2 _ n 

^( 0 ) — 1 ~ Z1 ai 1 0 ( 10 ) — 1 5 ( 10 ) — 

that play the role of obstructions to the embeddability of an abstract GL(2,M) structure into 
the Grassmannian Gr(3,5). It would be interesting to find a complete set of such obstructions. 

6 Concluding remarks 

This paper is a first step towards the general theory of integrability in Grassmann geometries. 
We gave a detailed characterisation of integrable systems S(X) associated with fourfolds X C 
Gr(3,5). 
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• We believe that some of our results can be generalised as follows. 

(a) In the dimension d = 3, the parameter space of non-degenerate integrable systems E(X) 
associated with submanifolds of codimension n — 3 > 2 in Gr(3, to) is finite-dimensional. 
Submanifolds X corresponding to ‘generic’ integrable systems are not algebraic. 

(b) In higher dimensions d > 4, every non-degenerate integrable system £(A) associated 
with a submanifold of codimension n — d > 2 in Gr(d, to) is necessarily linearly degener¬ 
ate. Submanifolds X corresponding to linearly degenerate integrable systems are rational 
(generally, singular). 

At the moment we are not aware of any (n, d)-independent approaches to the results of 
this kind. 

• Some examples of integrable systems (1) discussed in this paper can be interpreted as 
Backlund transformations: on elimination of v, they lead to a second-order PDE for u, 
similarly, on elimination of u they lead to a second-order PDE for v. It would be interesting 
to obtain a classification of Backlund transformations. 

• Our definition of integrability, based on the existence of holonomic trisecant submanifolds, 
applies to any abstract GL(2,M) structure. It would be of interest to understand whether 
every integrable GL(2,R) structure is necessarily embeddable, that is, comes from a four¬ 
fold X C Gr(3, 5). More generally, one may ask for a criterion of embeddability of an 
abstract GL(2,R) structure (see Remark at the end of Sect. 5.4). 

• It is a true challenge to classify integrable systems that correspond to algebraic fourfolds 
X C Gr(3, 5). The homology class of any such X can be represented as aa + brj where a, b 
are nonnegative integers, and a, rj are the standard four-dimensional Schubert cycles: a 
corresponds to 3-dimensional subspaces containing a fixed 1-dimensional subspace, and r) 
corresponds to 3-dimensional subspaces that have 2-dimensional intersections with a fixed 
3-dimensional subspace. Which values of a and b are compatible with the requirement 
of integrability? The approach of [9, 50] allows one to characterise algebraic X (in some 
special homology classes) as integral manifolds of certain overdetermined exterior differ¬ 
ential systems. Provided such characterisation is found for every a and b, it would be 
straightforward to intersect this differential system with our integrability conditions. 
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